The extent of introgression between incipient Clarkia species is determined by temporal environmental variation and mating system.

Introgression is pervasive across the tree of life but varies across taxa, geography, and genomic regions. However, the factors modulating this variation and how they may be affected by global change are not well understood. Here, we used 200 genomes and a 15-y site-specific environmental dataset to investigate the effects of environmental variation and mating system divergence on the magnitude of introgression between a recently diverged outcrosser-selfer pair of annual plants in the genus Clarkia. These sister taxa diverged very recently and subsequently came into secondary sympatry where they form replicated contact zones. Consistent with observations of other outcrosser-selfer pairs, we found that introgression was asymmetric between taxa, with substantially more introgression from the selfer to the outcrosser. This asymmetry was caused by a bias in the direction of initial F1 hybrid formation and subsequent backcrossing. We also found extensive variation in the outcrosser's admixture proportion among contact zones, which was predicted nearly entirely by interannual variance in spring precipitation. Greater fluctuations in spring precipitation resulted in higher admixture proportions, likely mediated by the effects of spring precipitation on the expression of traits that determine premating reproductive isolation. Climate-driven hybridization dynamics may be particularly affected by global change, potentially reshaping species boundaries and adaptation to novel environments.

Reproductive strategies and their consequences for divergence, gene flow, and genetic diversity in three taxa of Clarkia.

Differences in reproductive strategies can have important implications for macro- and micro-evolutionary processes. We used a comparative approach through a population genetics lens to evaluate how three distinct reproductive strategies shape patterns of divergence among as well as gene flow and genetic diversity within three closely related taxa in the genus Clarkia. One taxon is a predominantly autonomous self-fertilizer and the other two taxa are predominantly outcrossing but vary in the primary pollinator they attract. In genotyping populations using genotyping-by-sequencing and comparing loci shared across taxa, our results suggest that differences in reproductive strategies in part promote evolutionary divergence among these closely related taxa. Contrary to expectations, we found that the selfing taxon had the highest levels of heterozygosity but a low rate of polymorphism. The high levels of fixed heterozygosity for a subset of loci suggests this pattern is driven by the presence of structural rearrangements in chromosomes common in other Clarkia taxa. In evaluating patterns within taxa, we found a complex interplay between reproductive strategy and geographic distribution. Differences in the mobility of primary pollinators did not translate to a difference in rates of genetic diversity and gene flow within taxa - a pattern likely due to one taxon having a patchier distribution and a less temporally and spatially reliable pollinator. Taken together, this work advances our understanding of the factors that shape gene flow and the distribution of genetic diversity within and among closely related taxa.

Rapid evolution during climate change: demographic and genetic constraints on adaptation to severe drought.

Populations often vary in their evolutionary responses to a shared environmental perturbation. A key hurdle in building more predictive models of rapid evolution is understanding this variation-why do some populations and traits evolve while others do not? We combined long-term demographic and environmental data, estimates of quantitative genetic variance components, a resurrection experiment and individual-based evolutionary simulations to gain mechanistic insights into contrasting evolutionary responses to a severe multi-year drought. We examined five traits in two populations of a native California plant, Clarkia xantiana, at three time points over 7 years. Earlier flowering phenology evolved in only one of the two populations, though both populations experienced similar drought severity and demographic declines and were estimated to have considerable additive genetic variance for flowering phenology. Pairing demographic and experimental data with evolutionary simulations suggested that while seed banks in both populations probably constrained evolutionary responses, a stronger seed bank in the non-evolving population resulted in evolutionary stasis. Gene flow through time via germ banks may be an important, underappreciated control on rapid evolution in response to extreme environmental perturbations.

Phylogenomic analysis does not support a classic but controversial hypothesis of progenitor-derivative origins for the serpentine endemic Clarkia franciscana.

Budding speciation involves isolation of marginal populations at the periphery of a species range and is thought to be a prominent mode of speciation in organisms with low dispersal and/or strong local adaptation among populations. Budding speciation is typically evidenced by abutting, asymmetric ranges of ecologically divergent sister species and low genetic diversity in putative budded species. Yet these indirect patterns may be unreliable, instead caused by postspeciation processes such as range or demographic shifts. Nested phylogenetic relationships provide the most conclusive evidence of budding speciation. A putative case of budding speciation in the serpentine endemic Clarkia franciscana and two closely related widespread congeners was studied by Harlan Lewis, Peter Raven, Leslie Gottlieb, and others over a 20-year period, yet the origin of C. franciscana remains controversial. Here, we reinvestigate this system with phylogenomic analyses to determine whether C. franciscana is a recently derived budded species, phylogenetically nested within one of the other two putative progenitor species. In contrast to the hypothesized pattern of relatedness among the three Clarkia species, we find no evidence for recent budding speciation. Instead, the data suggest the three species diverged simultaneously. We urge caution in using contemporary range patterns to infer geographic modes of speciation.

Ancient Gene Duplications, Rather Than Polyploidization, Facilitate Diversification of Petal Pigmentation Patterns in Clarkia gracilis (Onagraceae).

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. Clarkia gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, that is, before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.

The opportunity for outcrossing varies across the geographic range of the primarily selfing Clarkia xantiana ssp. parviflora.

PREMISE: The timing of self-fertilization has potentially important consequences for the trajectory of mating system evolution, the opportunity for outcrossing, and the maintenance of genetic variation in populations. For primarily selfing taxa, it remains poorly understood as to how floral variation influences the opportunity for outcrossing and whether those attributes vary among populations across geographic ranges. METHODS: We examined variation in floral traits (herkogamy, protandry, flower size, stigma stage at anthesis, timing of stigma receptivity) in seven populations of Clarkia xantiana ssp. parviflora, a primarily selfing taxon, spanning from the western to eastern margins of its distribution. We also performed experimental emasculations and pollinations (followed by stigma severing) to quantify the extent of opportunities for outcrossing across flower development. RESULTS: There was marked among-population variation in all floral traits, particularly between far eastern and western populations. Emasculation experiments showed that the eastern populations had minimal autonomous selfing, but western populations had high rates of selfing within 24 h after anthesis. Population variation in autofertility was significantly predicted by floral trait variation, especially protandry and petal size. CONCLUSIONS: Greater protandry both extended the time over which outcrossing could potentially occur and reduced the probability of autonomous selfing, suggesting that there may be a tradeoff that results in fitness loss when pollinator visits are not common. The east-west pattern of differentiation in some floral traits parallels that of postglacial range expansion, suggesting that selection on the mating system may have been strong in the process of range expansion.

Floral density and co-occurring congeners alter patterns of selection in annual plant communities.

Although the evolution and diversification of flowers is often attributed to pollinator-mediated selection, interactions between co-occurring plant species can alter patterns of selection mediated by pollinators and other agents. The extent to which both floral density and congeneric species richness affect patterns of net and pollinator-mediated selection on multiple co-occurring species in a community is unknown and is likely to depend on whether co-occurring plants experience competition or facilitation for reproduction. We conducted an observational study of selection on four species of Clarkia (Onagraceae) and tested for pollinator-mediated selection on two Clarkia species in communities differing in congeneric species richness and local floral density. When selection varied with community context, selection was generally stronger in communities with fewer species, where local conspecific floral density was higher, and where local heterospecific floral density was lower. These patterns suggest that intraspecific competition at high densities and interspecific competition at low densities may affect the evolution of floral traits. However, selection on floral traits was not pollinator mediated in Clarkia cylindrica or Clarkia xantiana, despite variation in pollinator visitation and the extent of pollen limitation across communities for C. cylindrica. As such, interactions between co-occurring species may alter patterns of selection mediated by abiotic agents of selection.

Two genetic changes in cis-regulatory elements caused evolution of petal spot position in Clarkia.

A major premise in evolutionary developmental biology is that regulatory changes, often involving cis-regulatory elements, are responsible for much morphological evolution. This premise is supported by recent investigations of animal development, but information is just beginning to accumulate regarding whether it also applies to the evolution of plant morphology1-4. Here, we identify the genetic differences between species in the genus Clarkia that are responsible for evolutionary change in an ecologically important element of floral colour patterns: spot position. The evolutionary shift in spot position was due to two simple genetic changes that resulted in the appearance of a transcription factor binding site mutation in the R2R3 Myb gene that changes spot formation. These genetic changes caused R2R3 Myb to be activated by a different transcription factor that is expressed in a different position in the petal. These results suggest that the regulatory rewiring paradigm is as applicable to plants as it is to animals, and support the hypothesis that cis-regulatory changes may often play a role in plant morphological evolution.

Mating System Evolution under Strong Pollen Limitation: Evidence of Disruptive Selection through Male and Female Fitness in Clarkia xantiana.

Selection on floral traits in hermaphroditic plants is determined by both male and female reproductive success. However, predictions regarding floral trait and mating system evolution are often based solely on female fitness. Selection via male fitness has the potential to affect the outcomes of floral evolution. In this study, we used paternity analysis to assess individual selfing rates and selection on floral traits via male and female fitness in an experimental population of Clarkia xantiana where pollen limitation of seed set was strong. We detected selection through both female and male fitness with reinforcing or noninterfering patterns of selection through the two sex functions. For female fitness, selection favored reduced herkogamy and protandry, traits that promote increased autonomous selfing. For male fitness, selection on petal area was disruptive, with higher trait values conferring greater pollinator attraction and outcross siring success and smaller trait values leading to higher selfed siring success. Combining both female and male fitness, selection on petal area and protandry was disruptive because intermediate phenotypes were less successful as both males and females. Finally, functional relationships among male and female fertility components indicated that selfing resulted in seed discounting and pollen discounting. Under these functional relationships, the evolutionarily stable selfing rate can be intermediate or predominantly selfing or outcrossing, depending on the segregating load of deleterious mutations.

How petals change their spots: cis-regulatory re-wiring in Clarkia (Onagraceae).

A long-standing question in evolutionary developmental biology is how new traits evolve. Although most floral pigmentation studies have focused on how pigment intensity and composition diversify, few, if any, have explored how a pattern element can shift position. In the present study, we examine the genetic changes underlying shifts in the position of petal spots in Clarkia. Comparative transcriptome analyses were used to identify potential candidate genes responsible for spot formation. Co-segregation analyses in F2 individuals segregating for different spot positions, quantitative PCR, and pyrosequencing, were used to confirm the role of the candidate gene in determining spot position. Transient expression assays were used to identify the expression domain of different alleles. An R2R3Myb transcription factor (CgMyb1) activated spot formation, and different alleles of CgMyb1 were expressed in different domains, leading to spot formation in different petal locations. Reporter assays revealed that promoters from different alleles determine different locations of expression. The evolutionary shift in spot position is due to one or more cis-regulatory changes in the promoter of CgMyb1, indicating that shifts in pattern element position can be caused by changes in a single gene, and that cis-regulatory rewiring can be used to alter the relative position of an existing character.

Geographic variation in climate as a proxy for climate change: Forecasting evolutionary trajectories from species differentiation and genetic correlations.

PREMISE OF THE STUDY: Climate change models for California predict a warmer, drier future, potentially resulting in shorter growing seasons. If phenotypic differences between closely related species currently distributed across a moisture and temperature gradient represent adaptations to their abiotic environment, then as conditions become warmer and drier, populations presently adapted to cooler and wetter conditions may evolve to become more similar to those adapted to warmer and drier conditions. Two sister species, Clarkia unguiculata and C. exilis, are distributed across a moisture and temperature gradient in the southern Sierra Nevada, providing an opportunity to predict how this process may occur. METHODS: In a greenhouse experiment using wild-collected seeds from 11 populations in the southern Sierra Nevada, we examined relationships among elevation, climatic conditions, and population means for each trait, then evaluated bivariate relationships among maternal family means, using raw values and controlling for population and seed mass effects on phenotype. KEY RESULTS: Clarkia exilis occupied warmer, drier conditions, typically at lower elevations, than C. unguiculata did and flowered earlier and faster, producing smaller flowers with lower herkogamy. In C. unguiculata, petal area, herkogamy, and the rate of flower production were positively correlated with days to first flower. CONCLUSIONS: If selection favors earlier flowering, smaller petals, or faster flower production in C. unguiculata, then the genetic correlations among these traits should reinforce their joint evolution. Moreover, the correlations between these traits and herkogamy may promote the evolution of self-fertilization as an indirect response to selection, a previously unrecognized potential outcome of climate change.

Mating system divergence affects the distribution of sequence diversity within and among populations of recently diverged subspecies of Clarkia xantiana (Onagraceae).

PREMISE OF THE STUDY: The population biology of outcrossing and self-fertilizing taxa is thought to differ because of the advantage that selfers have in colonizing unoccupied sites where mates and pollinators may be limiting (Baker's Law). This reduced tendency for outcrossers to colonize new sites, along with their greater dependence on pollinators to disperse pollen, has the potential to differently influence the genetic diversity and structure of outcrossing and selfing populations. METHODS: We conducted a comparative population genetic study of two sister outcrossing and selfing subspecies of Clarkia xantiana that have very recently diverged. We used DNA sequence variation (>40 kb from eight nuclear loci) from large samples of individuals from 14 populations to assess geographic patterns of genetic diversity and make inferences about the demographic and colonization histories of each subspecies. KEY RESULTS: We show that sequence variation is strongly reduced across all selfing populations. The demographic history of selfing populations exhibits recent colonization bottlenecks, whereas such bottlenecks are rarely observed for the outcrosser. The greater effect of genetic drift in the selfer has resulted in strong population genetic structure, but with no pattern of isolation by distance. By contrast, the stronger effect of gene flow in the outcrosser has resulted in considerably less structure, but a significant pattern of isolation by distance. CONCLUSIONS: Taken together, our results suggest that selfing populations are not at migration-drift equilibrium, are affected by strong episodes of genetic drift during colonization, and experience little or no subsequent gene flow from other populations after those founder events.

Accelerated degradation of lignin by lignin peroxidase isozyme H8 (LiPH8) from Phanerochaete chrysosporium with engineered 4-O-methyltransferase from Clarkia breweri.

Free-hydroxyl phenolic units can decrease or even abort the catalytic activity of lignin peroxidase H8 during oxidation of veratryl alcohol and model lignin dimers, resulting in slow and inefficient lignin degradation. In this study we applied engineered 4-O-methyltransferase from Clarkia breweri to detoxify the inhibiting free-hydroxyl phenolic groups by converting them to methylated phenolic groups. The multistep, enzyme-catalyzed process that combines 4-O-methyltransferase and lignin peroxidase H8 suggested in this work can increase the efficiency of lignin-degradation. This study also suggests approaching the field of multi-enzyme in vitro systems to improve the understanding and development of plant biomass in biorefinery operations.

Rapid evolution of reproductive isolation between incipient outcrossing and selfing Clarkia species.

A major goal of speciation research is to understand the processes involved in the earliest stages of the evolution of reproductive isolation (RI). One important challenge has been to identify systems where lineages have very recently diverged and opportunities for hybridization are present. We conducted a comprehensive examination of the components of RI across the life cycle of two subspecies of Clarkia xantiana, which diverged recently (ca. 65,000 bp). One subspecies is primarily outcrossing, but self-compatible, whereas the other is primarily selfing. The subspecies co-occur in a zone of sympatry but hybrids are rarely observed. Premating barriers resulted in nearly complete isolation in both subspecies with flowering time and pollinator preference (for the outcrosser over the selfer) as the strongest barriers. We found that the outcrosser had consistently more competitive pollen, facilitating hybridization in one direction, but no evidence for pollen-pistil interactions as an isolating barrier. Surprisingly, postzygotic isolation was detected at the stage of hybrid seed development, but in no subsequent life stages. This crossing barrier was asymmetric with crosses from the selfer to outcrosser most frequently failing. Collectively, the results provide evidence for rapid evolution of multiple premating and postzygotic barriers despite a very recent divergence time.

Enhanced levels of S-linalool by metabolic engineering of the terpenoid pathway in spike lavender leaves.

Transgenic Lavandula latifolia plants overexpressing the linalool synthase (LIS) gene from Clarkia breweri, encoding the LIS enzyme that catalyzes the synthesis of linalool were generated. Most of these plants increased significantly their linalool content as compared to controls, especially in the youngest leaves, where a linalool increase up to a 1000% was observed. The phenotype of increased linalool content observed in young leaves was maintained in those T1 progenies that inherit the LIS transgene, although this phenotype was less evident in the flower essential oil. Cross-pollination of transgenic spike lavender plants allowed the generation of double transgenic plants containing the DXS (1-deoxy-d-xylulose-5-P synthase), coding for the first enzyme of the methyl-d-erythritol-4-phosphate pathway, and LIS genes. Both essential oil yield and linalool content in double DXS-LIS transgenic plants were lower than that of their parentals, which could be due to co-suppression effects linked to the structures of the constructs used.

Resource reallocation does not influence estimates of pollen limitation or reproductive assurance in Clarkia xantiana subsp. parviflora (Onagraceae).

PREMISE OF STUDY: Studies of pollen limitation and the reproductive assurance value of selfing are important for examining the process of floral and mating system evolution in flowering plants. Recent meta-analyses have shown that common methods for measuring pollen limitation may often lead to biased estimates. Specifically, experiments involving single- or few-flower manipulations per plant tend to overestimate pollen limitation compared to those involving manipulations on most or all flowers per plant. Little previous work has explicitly tested for reallocation within individual systems using alternative methods and response variables. • METHODS: We performed single-flower and whole-plant pollen supplementation and emasculation of flowers of Clarkia xantiana subsp. parviflora to estimate pollen limitation (PL) and reproductive assurance (RA). We compared levels of PL and RA using the following response variables: fruit set, seeds/flower, and seeds/plant. We also assessed the germination and viability of seeds to evaluate potential variation in pollen quality among treatments. • KEY RESULTS: Autonomous selfing in Clarkia xantiana subsp. parviflora eliminates pollen limitation and provides reproductive assurance. Estimates from single-flower manipulations were not biased, closely resembling those from whole-plant manipulations. All three response variables followed the same pattern, but treatments were only significantly different for seeds/flower. Pollen quality, as indicated by seed viability, did not differ among treatments. • CONCLUSIONS: Partial plant manipulations provided reliable estimates of pollen limitation and reproductive assurance. These estimates were also unaffected by accounting for pollen quality. Although whole plant manipulations are desirable, this experiment demonstrates that in some systems partial plant manipulations can be used in studies where whole-plant manipulations are not feasible.

Phylogeography of speciation: allopatric divergence and secondary contact between outcrossing and selfing Clarkia.

The origins of hybrid zones between parapatric taxa have been of particular interest for understanding the evolution of reproductive isolation and the geographic context of species divergence. One challenge has been to distinguish between allopatric divergence (followed by secondary contact) versus primary intergradation (parapatric speciation) as alternative divergence histories. Here, we use complementary phylogeographic and population genetic analyses to investigate the recent divergence of two subspecies of Clarkia xantiana and the formation of a hybrid zone within the narrow region of sympatry. We tested alternative phylogeographic models of divergence using approximate Bayesian computation (ABC) and found strong support for a secondary contact model and little support for a model allowing for gene flow throughout the divergence process (i.e. primary intergradation). Two independent methods for inferring the ancestral geography of each subspecies, one based on probabilistic character state reconstructions and the other on palaeo-distribution modelling, also support a model of divergence in allopatry and range expansion leading to secondary contact. The membership of individuals to genetic clusters suggests geographic substructure within each taxon where allopatric and sympatric samples are primarily found in separate clusters. We also observed coincidence and concordance of genetic clines across three types of molecular markers, which suggests that there is a strong barrier to gene flow. Taken together, our results provide evidence for allopatric divergence followed by range expansion leading to secondary contact. The location of refugial populations and the directionality of range expansion are consistent with expectations based on climate change since the last glacial maximum. Our approach also illustrates the utility of combining phylogeographic hypothesis testing with species distribution modelling and fine-scale population genetic analyses for inferring the geography of the divergence process.

Tempo and mode of mating system evolution between incipient Clarkia species.

Mating systems are among the most labile characteristics of flowering plants, with transitions frequently occurring among populations or in association with speciation. The frequency of mating system shifts has made it difficult to reconstruct historical evolutionary dynamics unless transitions have been very recent. Here, we examine molecular and phenotypic variation to determine the polarity, timescale, and causes of a transition between outcrossing and self-fertilization in sister subspecies of Clarkia xantiana. Phylogenetic analyses and coalescent-based estimates of the time to most recent common ancestor indicated that outcrossing is ancestral to selfing and that there has been a single origin of selfing. Estimates of divergence time between outcrossing and selfing subspecies were 10,000 (95% CI [credible interval]: 3169-66,889) and 65,000 years ago (95% CI: 33,035-151,448) based on two different methods, suggesting a recent and rapid evolutionary transition. Population genetic data indicated that the transition to selfing was associated with a 80% reduction in molecular diversity, which is much greater than the 50% reduction expected under a shift from obligate outcrossing to obligate self-fertilization alone. Our data also suggest that this severe loss of diversity was caused by colonization bottlenecks. Together with previous studies, evidence for reproductive assurance in C. xantiana now connects variation in plant-pollinator interactions in the field to phenotypic and molecular evolution.

Physiological performance and mating system in Clarkia (Onagraceae): does phenotypic selection predict divergence between sister species?

PREMISE OF THE STUDY: The evolution of self-fertilization often occurs in association with other floral, life history, and fitness-related traits. A previous study found that field populations of Clarkia exilis (a predominantly autogamous selfer) and its sister species, Clarkia unguiculata (a facultative outcrosser) differ in mean photosynthetic rates and instantaneous water use efficiency (WUE(i)). Here, we investigate the strength and direction of selection on these traits in multiple populations of each taxon to determine whether natural selection may contribute to the phenotypic differences between them. METHODS: In spring 2008, we measured instantaneous gas exchange rates in nine populations during vegetative growth (Early) and/or during flowering (Late). We conducted selection gradient analyses and estimated selection differentials within populations and across pooled conspecific populations to evaluate the strength, direction, and consistency of selection on each trait early and late in the season. KEY RESULTS: The direction and relative strength of selection on photosynthetic rates in these taxa corresponds to the phenotypic difference between them; C. exilis has higher photosynthetic rates than C. unguiculata, as well as stronger, more consistent selection favoring rapid photosynthesis throughout the growing season. Patterns of selection on transpiration, WUE(i), and the timing of flowering progression are less consistent with phenotypic differences (or lack thereof) between taxa. CONCLUSIONS: We detected several examples where selection was consistent with the phenotypic divergence between sister taxa, but there were also numerous instances that were equivocal or in which selection did not predict the realized phenotypic difference between taxa.

Population genetics and the evolution of geographic range limits in an annual plant.

Abstract Theoretical models of species' geographic range limits have identified both demographic and evolutionary mechanisms that prevent range expansion. Stable range limits have been paradoxical for evolutionary biologists because they represent locations where populations chronically fail to respond to selection. Distinguishing among the proposed causes of species' range limits requires insight into both current and historical population dynamics. The tools of molecular population genetics provide a window into the stability of range limits, historical demography, and rates of gene flow. Here we evaluate alternative range limit models using a multilocus data set based on DNA sequences and microsatellites along with field demographic data from the annual plant Clarkia xantiana ssp. xantiana. Our data suggest that central and peripheral populations have very large historical and current effective population sizes and that there is little evidence for population size changes or bottlenecks associated with colonization in peripheral populations. Whereas range limit populations appear to have been stable, central populations exhibit a signature of population expansion and have contributed asymmetrically to the genetic diversity of peripheral populations via migration. Overall, our results discount strictly demographic models of range limits and more strongly support evolutionary genetic models of range limits, where adaptation is prevented by a lack of genetic variation or maladaptive gene flow.